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Abstract 

Using a three-baryon ANN potential having a spin-independent dispersive 
and a two-pion exchange components, realistic variational Monte Carlo cal- 
culations have been performed for the A-seperation energy for ^He. A rela- 
tionship between strengths associated with both pieces of ANN force giving 
exact experimental A-seperation energy has been obtained which then com- 
bined with our previous study on hypernucleus determines a unique and 
acceptable set of strength parameters for ANN potential. The ANN force is 

found important for the core polarization. 

PACS numbers: 21.80.+a,21.10.Dr,13.75.Ev,27.20.+n 
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The Cluster Monte Carlo(CMC) technique first used for nuclei |l] is recently generalised 
to hypernucleus for a realistic microscopic calculations @. The same has been further 
extended to A He in this work to study the A-dependence of the expectation value of ANN 
potential which is of great importance. 

In case of ^O, we find that the operatorial NN, NNN, AN and ANN correlations 
completely change the previously determined ANN potential using central correlations |§ 
and that a reasonable Vajvtv results in a A-seperation energy(I?A) consistent with the "em- 
pirically" determined value i.e. -13 ±0.4 MeV. Also, ANN force generates a corresponding 
correlation which has been found playing a significant role in changing the density profile 
and energy of the core nucleus. This provokes a deep interest in determining the correct 
three-baryon AiViVpotential. 



The ANN potential is written as a combination of a dispersive central force [[3] (V/° 



ANN) 



r2n \ 
'ANN) 



and a two-pion exchange force f|] (V^ 

Vann = Vann + v ann- (!) 

We do not invoke a spin dependence in the V ANN which is unimportant for spin-zero core 
hypernuclei. Thus, we use a phenomenological form 

VSnk = W T*(r lA )T*(r 2A ). (2) 

Indeed, we find a negligibly small value for spin-spin AA" potential in and also herein 
for ^He. T n (r) is the OPE tensor potential 

T.= fl + - + i)-(l-e-f (3) 



with x=/ir, fi=0.7 fm -1 is the pion mass, and the cut-off parameter c=2.0 fm -2 . The V A ^ fN , 
consists of two parts, arising from s- and p- wave 7r — A potentials Q . The former is strongly 
suppressed and that almost all contribution is coming from the latter one. We, thus, do not 
consider the s-wave channel and write the V A ^ N as 

V ANN = ~ ( r i ' T 2) {ci ■ cta^tt^ia) + S 1A T n (r 1A ), cr2.(T A Y n (r2 A ) + S 2 AT 7r (r 2A )} , (4) 



The SiA is the tensor operator and 

YJr) 



-/j/r 



fir 



(5) 



The experimental B A value has been obtained theoretically by the relation 

< V N \H N \V N > _ < jjijjjr > 
A < V N \V N > < V\V > ' U 

where, \I/ is the full wavefunction of the hypernucleus, \I/jv is the ground state wavefunction 
of A-l nucleons. The H and Hn are the hamiltonians for the hypernucleus and its core, 

-2 A-l A-l 



H = H N — 



Tr 



2m\ 



(7) 



i=l 



A-l fi 2 A-l A-l 

^ = -E|rV? + E% + E *W 

i=l i<j i<j<k 



(8) 



where and Vijk refer to NN Argonne potential H and NNN Urbana model VII potential 
||. The two-baryon Van force is written as a sum of purely central, space exchange and 
spin-spin pieces as 

1 



^Aiv(r) = u (r)(l - e + eP a .) + ^; CT T^(r)cr A ■ a N , 



(9) 



vo(r) = v c (r) - vT%(r), 



(10) 



where P x is the space exchange operator and e is the corresponding exchange parameter. 
The v = (v s + 3u t )/4 and t> CT = t> s — v t are respectively the spin-average and spin- dependent 
strengths, where v s and v t denote the singlet and triplet state depths, respectively Fi- 
nally, f c (r) is a short range Saxon-Wood repulsive potential. The various parameters are 
w s =6.33MeV, v t =6.1MeV, e =0.3, W c =2137MeV, R=0.5fm, anda=0.2fm. The varia- 
tional wavefunction for |He under the independent triplet approximation is written as 



n (i + u iA ) 

Li=l 



sl[(i + u i: 

i<3 



i + E^-a+ E u ijk 

i<j i<j<k 



3 



A-l 

n 

i=i 



A-l 



i*j >= n / c a (^a) n /c^oai^- 1 > 



(12) 



The product A^" 4 1 > is a Slater determinant, where A is the antisymmetrization operator 
running over nucleons 



Al**" 1 >= A [(n T)i(n |) 2 (p t) 3 (p i) 4 (A |) 5 ] 



(13) 



The three-baryon correlations 



(14) 



have the same structure as the three-baryon potentials though it differ through the range c 
of the cut off functions T w (r) and Y w (t). The subscript B stands for N and A both. For 
and Ui\ we make the choice 



Uij — f3pUp(rjj)Ofj , 

p=2 



(15) 



C/ iA = a CT w A (r; A )cr A • <r h 



(16) 



with 



f A f A 

A _ if is 

~~ f A 



/i 



The / A is the spin averaged correlation function 



(17) 



(18) 



and / A and / t A are the solutions of Schrodinger equations with quenched AN potentials in 
singlet and triplet states respectively 

h 2 



V +v s(t) (r AN ) + 9 AN 



fit) = o. 



(19) 



The quenched potential v styt ) is quenched in the two-pion parts of the central and spin 
channels 
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v s (r) = a a v c (r) - (a 27r v + ~a a v a )T*(r), (20) 



v t (r) = a a v c {r) - (a 2n v - ^a a v a T%(r), (21) 



and 



Ti 2 ( 2 2k an (u an - 1) v kN (v kN - 1)\ A ^ 

#A7V = " I «A7V : + ~ 2 1 - e 



2//AJV V 

7a;v 
1 + e a Aiv 



2 • 
AN 



The ^AAf is such that the A-nucleon correlations in eq(ll) has the asymptotic behaviour 
required by the A-l body Schrodinger equation 
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f C AN ~ r—e~^ r . (23) 

Also in Uij, a central three-body correlation is folded as 

Ua - II kkUiy (24) 

We refer to ref. @ for details of fijk, Unn and Unnn correlation functions. 

First of all, we made variational searches for 4 He of ref. and find that every parameter 
turns out to be the same except the which is -.005 instead of -.004. Then we made 
variational parameters searches for A He for many combinations of Wq and C p including 
Wo = C p = (which means Uann an d Vaatat have been set off), we find that a a is unity 
, however, unlike ^O, a 2n is found sensitive to C p obeying 

a 2n = 0.92 -0.1 C p , <C P < 0.7 MeV. (25) 

We also note an increase in a 2n from its minimum value at C p = 0.7 

a 2n = 0.85 + 0.2 C p , 0.7 < C p < 1.0 MeV . (26) 

This is probably because U ann and the repulsive / A both go together to suppress the 
nucleon density in the interior region unless the most stable position is achieved and then 



f A gets increased to cancel the effects of Uann with further enhancement of C p . In case 
of this sentivity was weak. The parameter ea is found to be -.0014 and the cut off for 
Uann is 1.6 fm -2 . The variational parameters used in eq.(22) are KAJV=0.02fm, van=0-7 
Can =2.3 fm~ 2 , aAAr=1.4fm, rAAr=0.7fm, and jan is eigen value parameter determined 
by matching logarithmic derivatives at a suitable r. The 4 He energy in our calculation 
corresponds to the first six operators of Argonne Vu potential and is 31.1±0.1 MeV. The 
CMC technique up to the four-body clusters is complete for the core with four nucleons 
,however, we neglect a sole five-body cluster involving a A. 

We observe that V A 2 A y N is substantial and negative and that it is due to the presence of 
non central correlations in the wave function as its values are -2.65 ±0.05 and 0.19 ± 0.03 
with and without non central correlations, respectively. We also note that S\% operator in 
U\2 changes the < V^f 2 >• These results are along the line of ^O. 

We find a strong dependence of Ba values on the ANN potential which has been fit in 
a functional form depending upon the strengths Wq and C p as 

B A (iZ) =a + bC p + cC 2 p + dW . (27) 

The dependence on Wq is linear, however, the quadratic dependence on C p is attributed 
to the Uann correlation. The constants a, b, c and d depend upon the hypernucleus j^Z. 
The constant a refers to the theoretically calculated -Ba value when V 2 ^ N force and its 
corresponding correlation is set off. We variationally determine Ba for many combinations 
of W and C p . For example, for Wq = C p = 0, B A (%H.e) is 4.8 ± 0.1 and for Wq = .016, C p = 
0.7, B a is 2.9 ± 0.1. we then perform correlated difference Monte Carlo calculations using 
the Wo = .016, C p = 0.7 random walk for other sets of Wq and C v . We fit all these Ba and 
5Ba (given in table I) values using above functional form and get a=4.8 MeV, b=2.45, c=2.5 
MeV -1 and d=-291. Using these constants, one may then obtain a relationship between Wq 
and C p for the experimental -Ba (|He) i.e. 3.12±0.02 MeV ||. A similar relationship between 
Wq and C p has already been obtained for hypernucleus for its empirically etermined -Ba 
value. We then plot Wq and C p for both the hypernuclei and get a unique combination of Wq 
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and C p where the two plots intersect each other (fig. 1). This determines the ANN potential 
with strengths Wq = 0.015 MeV and C p = 0.67 MeV. Using these strengths, we reperform 
the calculations and get £?AQHe)=3.1±0.1 MeV. The core energy is -23.4 ± 0.2 MeV in this 
case. However, it is -21.0 ± 0.2 MeV with no U\nn and V\nn in the calculation. Thus a 
distortion of 2.4 MeV is found due to ANN potential. The results are given in Table II and in 
Table III for the newly determined potential and for Wq = C p = 0, respectively. Accidently, 
we chose Wo = -015 and C p = 0.7 for the in our previous study which are very close 
to the newly determined strengths and had got the _Ba(a 7 0)=13.5±1.8 MeV slightly higher 
than the empirical B^. Certainly, a decrease in C p from 0.7 to 0.67 will reduce the negative 
contribution of Vl^ N and will bring the results very close to the emperical value. Thus, the 
conclusions drawn over there are correct for the new Vaaa potential too. 

It is obvious from fig. 2 (where point nucleon and A densities are shown) that A is strongly 
localised in the interior region of the core and that it pushes away the nucleons towards lower 
density regions both periphery and centre of the nucleus, presumably due to repulsive /ajv 
correlation. The £/aaa reduces the A density in the interior region and thus as a result 
nucleons are less pushed away. In both cases of nucleon and A densities, Vaaa and its 
corresponding correlation Uann are found to play a significant role. 

The addition of the space-exchange correlation in L^a (eq.16) is very unlikely to change 
the situation as it is weak and that the cut of the curves in fig.l is sharp. We are in a way to 
include this to see its effect specially on the hypernucleus where many p-wave nucleons 
are present. Also, the use of spin dependent V® NN for spin non-zero core hypernuclei will 
be interseting. Such a study for ^Li, A He and is in progress. We will also try to include 
the s-wave part of the V\nn- Moreover, we believe that U\nn should be responsible for the 
A spin-orbit splittings in various excited states of hypernuclei as there is no two-body AN 
tensor correlation. Using the newly determined Cp we will be able to calculate it for the 
for the A being in p-shell with J=3/2 and 1/2. The experimental B\( l ^O) for A being 
in p-shell is 2.5 ± 0.5 MeV [[|. Thus the splitting is less than 1 MeV. In other words, this 
study will provide a check on C p . 
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FIGURES 
FIG. 1. C p vs W . 

FIG. 2. One-body densities for the nucleons and A in ^He, and for 4 He. 
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TABLES 

TABLE I. Differences B A (C P , W )-B A (C P = 0.7, W = 0.016). 



C p 


Wo 


5B A MeV 


0.7 


0.02 


-0.98 ± 0.03 


0.7 


0.018 


-0.49 ±0.01 


0.7 


0.014 


+0.5 ±0.02 


0.8 


0.02 


-0.46 ± 0.07 


0.8 


0.016 


±0.71 ±0.05 


0.8 


0.014 


±1.31 ±0.05 


1.0 


0.022 


±0.21 ±0.1 


1.0 


0.018 


±1.54 ±0.1 


1.0 


0.016 


±2.22 ±0.1 
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TABLE II. Variational Energies for C p = 0.67 MeV, W = 0.015 MeV. All values are in MeV. 
Numbers in parentheses are statistical errors in the last digit. 



Clusters 




one- body 


two- body 


three- body 


iour-body 


lotal 


A Kinetic Energy a 




11.6(1) 


0.02(1) 


1.08(4) 


0.03(2) 


12.76(8) 


A AT T~l j. j • 1 

AJN Potential 


v (r)(l - e) 




-17.10(9) 


( ~\ fin ( c\\ 

-0.63(2) 


-0.03(3) 


1 f~7 'If* ( 1 f\\ 

-17.76(10) 




v (r)eP x 




'I f\f\ ( A \ 

-7.09(4) 


0.50(1) 


0.09(1) 


-6.50(4) 




1 n-<2 




0.09 


-0.03 


0.01 


0.07 


AJNJN Potential 


'ANN 






3.23(3) 


0.11(1) 


O O A ( 0\ 

3.34(3) 




T/27T 

ANN 






-2.30(3) 


-0.36(2) 


-2.66(4) 


A TV 

A Energy 




11.6(1) 


-24.1(2) 


i n/i \ 

1.9(1) 


-.15(3) 


-10.7(1) 


Nuclear Kinetic 




81.1(3) 


47.0(1) 


7.6(1) 


-.3(1) 


135.4(4) 


NN Potential 


Vij 




-144.9(3) 


-3.5(1) 


1.2(1) 


-147.2(3) 


NNN Potential 


Vijk 






-11.3(1) 


-.28(3) 


-11.6(1) 


Nuclear Energy 




81.1(3) 


-97.9(2) 


-7.3(1) 


.6(1) 


-23.5(1) 


Total Energy 




92.7(3) 


-122.0(3) 


-5.4(1) 


.4(3) 


-34.2(1) 



a Includes nucleon kinetic energy from A correlations. 
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TABLE III. 


Variational energies for C p 


= W = 0. 


All values are 


in MeV. Numbers in 


parentheses are s 


statistical errors in the last d 


igit. 








Clusters 


one-body 


two-body 


three-body 


four-body 


Total 


Kinetic Energy a 


13.2(1) 


0.05(1) 


0.34(4) 


-0.09(2) 


13.5(1) 


AN Potential 


v (r)(l - e) 


-19.4(1) 


-0.72(4) 


0.03(5) 


-20.1(2) 




v (r)eP x 


-8.2(1) 


0.56(2) 


0.18(2) 


-7.5(1) 




\v a T%a A ■ <t n 


0.1 


-0.04 


0.01 


0.1 


A Energy 


13.2(1) 


-27.5(2) 


0.1(1) 


0.1(1) 


-14.0(1) 


Nuclear Kinetic 


82.9(4) 


49.5(2) 


8.0(2) 


-0.2(1) 


140.2(5) 


NN Potential 


Vij 


-147.5(4) 


-4.1(2) 


1.4(1) 


-150.1(4) 


NNN Potential 


Vijk 




-11.7(1) 


-3.2 


-12.0(1) 


Nuclear Energy 


82.9(4) 


-97.9(3) 


-7.8(1) 


0.8(1) 


-21.0(2) 


Total Energy 


96.1(4) 


-125.4(5) 


-7.7(1) 


0.9(1) 


-35.9(1) 



a Includes nucleon kinetic energy from AN correlations. 
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